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Abstract 

This paper is devoted to the convergence analysis of the Adini element 
scheme for the fourth order problem in any dimension. We showed that, 
under the regularity assumption that the exact solution is in H*, the 
Adini element scheme is 0{h?) order convergent in energy norm, and the 
convergence rate in norm can not be nontrivially higher than 0{h^) 
order. Numerical verifications are presented. 
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1 Introduction 


This paper is devoted to the convergence analysis of the Adini element 
scheme for the fourth order problem in any dimension. The Adini element, (c.f. 
[H for 2D, m for higher dimension) is among the earliest finite elements for 
elliptic problems. It uses the rectangles (2D) and generalised rectangles (higher 
dimensions) as geometry shapes, and the evaluation and the derivatives of first 
order on the vertices as nodal parameters. The generation of stiffness matrix is 
easy and friendly, and this element has become a popular one during the past 
half century, and stimulated various works, (see 0, m, M)- In this 

present paper, we discuss the capacity of the Adini element scheme for fourth 
order problems, and present a sharp analysis of the upper and lower bound of 
the convergence rate in energy and integral norms in arbitrary dimensions. 

When used for second order problems, the Adini element scheme is a con¬ 
forming one, and the error analysis is straightforward by the fundamental Cea 
lemma and standard arguments. When used for fourth order problems, however, 
the Adini element is a nonconforming one, and the convergence analysis is more 
subtle. In Wang, Shi and Xu [14] where they generalised the Adini element from 
2D to arbitrary dimension, the 0{h) convergence rate of Adini element in any 
dimension has been proved for fourth order problem. Meanwhile, a higher accu¬ 
racy of the scheme is still expected and numerically observed. In 2D, it has been 
proved by Lascaux and Lesaint [5] that the finite element solution converges to 
the exact solution with 0{h^) order in the energy norm, provided the rectan¬ 
gular cells in the grid are all the same. Then in 2004, Lin and Luo nni showed 
the 0{h?) convergence of the Adini element without assuming the congruence 
of the cells of the grid. Later in 2006, Mao and Chen m showed further the 
0{h'^) convergence rate for anisotropic grids. So far to our knowledge, the sharp 
analysis of the convergence rate of Adini element for fourth order problems in 
higher dimensions is still absent. 

In this paper, we study the convergence rate of the Adini element scheme 
for fourth order problems in higher dimensions. Technically, without making 
a crucial use of the nodal interpolation which was done by (0, uni, m) and 
which will bring extra regularity assumption on the exact solution in higher 
dimension, our analysis relies on the structure of Adini element space only. We 
figure out the intrinsic symmetry property of the Adini element space, and show 
the 0{h?) energy norm convergence rate in a unified way with respect to the 
dimensions provided the exact solution belongs to H'^. 

There have been works that study high accuracy nonconforming finite ele¬ 
ment methods for fourth order problems. Several 0{h?‘) nonconforming finite 
elements have been constructed in, e.g., 0(11 ng. In contrast to these ele¬ 
ments, the Adini element space does not possess such moment continuities; the 
average of the normal derivatives of Adini element function is not continuous 
across the internal faces. This hints us to make use of a different way by using 
the symmetric property inside one cell, and moreover, this unusual property 
makes it hardly possible to make use of the dual argument to obtain higher 
order convergence rate in or norm. Indeed, in the paper, we further 
show that the convergence rate in norm can not be non-trivially higher than 
0{h?) order. 

The analysis of the lower bound of the convergence rate of the Adini element 
scheme in norm is a generalization of Hu-Shi’s work 0, which solved an open 
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problem whether the convergent order in norm can always be higher than 
that in the energy norm. Technically, a decomposition of the residual (/, u — Uh) 
to a leading term and other higher order terms works crucially, and we estimate 
the lower bound of the leading term sufficiently. Again, a sharp analysis of the 
interpolation operator will play a key role. Therefore, by the discrete Poincare 
inequality, we obtain that the convergence rate of the Adini element scheme for 
fourth order problem in energy norm, norm and the integral norm are all of 
O^h?) order, and these estimates are all sharp. 

The remaining of the paper is organized as follows. In Section [21 we present 
some preliminaries of the Adini element. In Section l3.ll we present the model 
problem and the Adini finite element discretization. In Section [3.21 we show the 
0{h^) order convergence rate in energy norm in any dimension. In Section (3.31 
we further show the 0(h?) order convergence in norm in any dimension. In 
Sectional some numerical examples are presented to demonstrate our theoretical 
results. Finally, in Section 0 some conclusions are given. 

2 The preliminaries: the Adini element 

2.1 The Adini element 

Let K C be a d-rectangle, Xc = (a;i,c,a: 2 ,c, • • • ,Xd,c)^ G be the 
barycentre of AT, and hi the half length of K in Xi direction, i = l,2,...,d. 
Then the d-rectangle can be denoted by 

AT = {x = (xi, X 2 , • • • , Xd)^ I Xi = Xi^c + iihi, -1 < < 1, 1 < z < d}. 

Particularly, the vertices Oi, 1 < z < 2^^, of AT are denoted by 

X2,c~h^i2 ^2 , * ‘ , Xd^c~h^idhd') 7 I — ^ — j — 1 ^Z ^2 . 

Moreover, denote by and ^ the two {d — l)-dimensional faces of K 
without the edges parallel to the Xi axe; see Figure 1. 

The d-dimensional Adini element is dehned by the triple {K, Pa{K), D), 
where 

• the geometric shape AT is a d-rectangle; 

• the shape function space is 


Pa{K) := Qi{K) -k span{x-g | 1 < z < d, q€ Qi{K)}, (1) 


here and throughout this paper, Qi{K) denotes the space of all polynomi¬ 
als which are of degree< I with respect to each variable Xi, over AT; 

• the nodal parameters are, for any v € C'^(Ar), 



where are vertices of AT, z = I,..., 2^^. 


3 




Figure 1: degrees of freedom for the Adini element 


Let a denote the multiple-index with a = (ai, • • • ,ad), Oii{^ < i < d) are 

d d 

nonnegative integers, and |a| = ^ a^, x°‘ = Y[ ■ The partial derivative 

i=l 1=1 

operator can be written as 


9“ = 






Let Ci {1 < i < d) be the d—dimensional unit multi-index with its i—th entry 
equal to 1. 


2.2 Structural properties of the Adini element 

Given a d-dimensional domain, 7/i is a triangulation on tl, and K G 
Til- Let be piecewise bilinear interpolation operator on K, namely G 
Qi{K) and (n)^u)(P) = v{P), for any vertex P of K, and v G C(K). Define 
on C{K) the operator TZ\. := Id — II)^, with Id being the identity operator. 

Define Ilo.icw = wdx, for any w G lf{K). The global version IIo of the 

interpolation operator Hq^k is defined as 

Itolif = IIo^K, for any K gTh- (3) 

Lemma 2.1. It holds for wt G Pa{K) that 

, l<i<d. (4) 

Proof. Given Wh G Pa{K), a direct calculation leads to that 

€ Qi{K) + span {(xj - Xj^cf • d, d e Q\{K), 1 < j < d} , (5) 

where Q\{K) := span{(a; - Xc)°'}\a^\<i^ai=o- Further, 

S]i=s-paxi{{{xj-Xj^cf-h'^)-q,l<j<d,qGQ\{K)]. ( 6 ) 


Tpl dwh 


= n 


dwh 


K 


dxi 
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Noting that {xj — Xj^c)'^, I < j < d, evaluate the same on ■ and F'j^ ■, we 
obtain This finishes the proof. 

□ 


For ease of expression, we define the following sets 

Mi j — {(cKi,*'* ,o;c/)| crj — 1, 2 ^ ckj ^3, k ^ ^ 

= {(ai.a 2 , 0 : 3)1 ai = 1, 2 < < 3, as, < 1, 2 < k ^ j < 3}. 

By means of ([51), on F'^ ■ of the element K, we can get that 
,1 dwh' 


n 


K 


dxi 


(^ 1 ; ‘ ; ^i—1 7 ^i,c i ; ^z + 1 ; 


,Xd) 


= Bf{j,a)Uo,K{d°‘wh)7 


l<j<d a^Mij 

3¥=i 


with 


ddi (ji o) — ,\{xj Xj^c) ■’ dj{xj Xj^c) ^ ] (a; Xc) 


oc—ei—ct-j e-i 


ajl 


d°‘wh G span{l,Xj}, if aj = 2, and d'^Wh is constant, if aj = 3. 


Noticing that {j,a) = Q,iii^ j, and flo^ic (9“w/i) are constant. 


(7) 

( 8 ) 
(9) 


For example, in two-dimensional case, 

1 .. ^3 7 2 2 9'^Wh 

+ gl(7:7-777,J 

and in three-dimensional case 

F](^^{xi^c ± hi,X2,X3) 

3 1 

= ^ —[{xj - Xj,cT^ - h^jixj - Xj,cr^-^]{x - Wh) 

j=2a'eM[ . 

Given K G Fh, we define the canonical interpolation operator Hk '■ C^{K) —>• 
Pa{K) by, for any v G 


{Ukv){P) = v{P) and (VnKi;)(F) = Vn(P), 


( 10 ) 


for any vertex P of K. The interpolation operator 11^^ has the following error 
estimates: 

\v-IIkv\i,k / = 0,1,2,3,4, (11) 

provided that v G F[^{K), where s > 4 and s > | -|-1 such that F[^{K) C C^{K), 
see Remark 13.41 
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Lemma 2.2. For any u € Pa{K) and v G Pa{K), it holds that 


{v\u-nKu),v\)L^K) = -Y. E y 

i=l l<i<d ^ 

Proof. It follows from the definition of Pa{K) that 
d'^v 

d^v 


h? f ay d‘^1 


IK 


dx^dx] dx} 


dx, (12) 


€ Qi{K)y 


dxidxj 


G{p-q\pGPi{xi,Xj), qGQYiK)} 

+ span{xfc -q, qG Q^iK), 1 < fc < d}, i ^ j, 


(13) 


where Pi(xi,Xj) := spanjl, Xi, xj}, Qi\K) := span{x“}|a^|<i,a,=a^.=o- 
Since u G Pi{K), we have, with 


hf^d^u. “ ay ^ 2^2 

“=“■ + it E ^ E E 

1=1 * 1=1 l<j<d * 1 

j>i 

where ui G Pa{K)- 

The Taylor expansion and the definition of the operator 11^^ yield 


hi ’ 


d 


(14) 


i=i * i=i i<j<d * 1 

j>i 


d 


Thus 




axf 

a^(M — n^u) 


ay 


dxidxi 




Since 


and 


[ - 4:)qkdx = 0, gfe e Qi{K), 1 < k i < d, 

JK 


(16) 


[ (c| - l)^idx = 0, / (^1 - l)cyx = 0, / (^1 - l)C*?idx = 0, 1 < t 7^ j < d, 

JK JK JK 

(17) 

a combination of (USD and m and some elementary calculation yield 


d‘^{u — Ilii-u) d'^v 


3 .Ik dxfdxj dxf 


Jk dxf dxf 
By the same argument, it yields 

f d‘^{u — li-Ku) d'^v 


, h^ f d'^u d'^v , / nx 

dx = —y / „ n „ n dx, 1 < I 7 ^ j < d. (18) 


dxidxj dxidxj 

which completes the proof. 


dx = 0, 


(19) 

□ 
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3 The capacity of Adini element for fourth order 
problems 


3.1 Model problem and finite element discretisation 

Let 12 C be a bounded domain with Lipschitz boundary. In this paper, 
We consider the model fourth order elliptic problem: 


f = /, in 12 , 

\ M = 0, on an. 


( 20 ) 


The variational formulation is, given f G H ^(12), to find u € V := ilQ(r2), 
such that 

an{u,v) = (/,u), for any v € V. (21) 



Let Th he a regular d-rectangle triangulation of the domain 12. Define the 
Adini element space in a standard way by 

14 := {vGL'^in) ■ v\k e Pa{K), vat eTh, V and Vv is continuous at 

all internal vertices}, 

and associated with the boundary condition, 

14o := {vh G 14 : Vfi and Vvh vanishes at all boundary vertices}. 

Evidently, 14 C H^{Q) and 14o C ilo(12) ([13]). However, 14 and 

14o ^ Evidently, P^{K) C Pa{K) for any K G Tn- By the standard 

technique, 

inf |u-u^|z,^<C/i4-'|u|4.n, ; = 0,1,2,3,4, (22) 

vheVho 

for any v G iL^(12). Herein and throughout this paper, C denotes a generic 
positive constant which is independent of the meshsize and may be different at 
different places. 

Associated with the model problem, the Adini finite element problem is to 
find Uh G 14o, such that 

ah{uh,Vh) = {f,Vh)L'^(Q), for any Vh G 14o, (23) 

where ah{uh,Vh) ■■= ^ aK{uh,Vh). 

KeTh 

Define a semi-norm over 14 by \uh\\ '■= J^KeTh Poincare 

inequality, | • |/i is a norm on 14oj and it is equivalent to || • while the latter 
denotes the piecewise norm. 

3.2 Error analysis in energy norm 

In this section, we present an upper bound of the energy norm of the error 
of the hnite element scheme (123. A main result of this paper is the theorem 
below. 


7 


Theorem 3.1. Let u and Uh he the solutions of o and (1231); respectively. 
Assume that u G Then 

\u — Uh\2,h ^ Ch'^\u\A.Q.- (24) 

Proof. By the second Strang Lemma, we have 

\u-Uh\h<C{ inf \u-v\h+ sup ( 25 ) 

■«^Vho WhGVho I'^hlh 

where 

Eh{u,Wh) ■■= ah{u,Wh) - {f,Wh) = ^ (26) 

The first term of (1251) is the approximation error and the second one is the 
consistency error. 

We shall consider separately the faces orthogonal to the xt axes (1 < * < d), 
namely we rewrite the consistency error to 


Eh{u,Wh) = y^^E^.{u,Wh), 


(27) 


with 


E^.{u,wh) = 

K&Th 

= E 

KgTh 


f d'^u dwh 

IdK dxi 


n^ids 


’dK 


dn^ dx, 


KGTh 


dWh 
dxf ^ dxi 


ds 


d^i 


.- V v} 


K&Th 


dwh 
dxf ’' dxi 


(28) 


where is the unit outward normal parallel to the Xi axe. 








Let K ^Th and denote g = It holds that, by (O, 


^ K,i ^ K,i 


^K,i l<j<daGMij 


5 X] X Bf{j,a)Ilo,Kd°‘whds 

<3<^ 


[ ^ 

Ik dxi 

f ^ 

Ik dxi 


X X -Bf (j,Q:)no.K5“w/,dx 


l<j<d ct^Mij 


(29) 


X X Bf(j^O')d°‘whdx 


l<j<d a^Mij 

j¥^i 


+ / ^ X X B^i.j,oi){kVQ^K - Id)d°‘whdx 




:= Lfi +L5. 


Integrating by parts yields 

Ik 


1-'^ ^ ^ Sf (j,a)a“M;,ida; 


d'^g 


Ik dxi 


X X Bi^{j,a)d°^ ^'Whdx 


l<j<d a^Mij 


(30) 


+ 


iL S S Bf(i,o)8‘-.»id,^ 

■^^K.i'' l<j<daeMij 

Since u € iL^(I2), ic/j e H^{il) and d°‘~^*Wh, (a G IWij) are tangential deriva¬ 
tives of the faces that orthogonal to the axe Xi, thus and d°‘~'^'Wh, (a G Mij) 
are continuous across faces F'k^^ F'l,. we obtain 




X ^M = - X / 7^ X X FfU,a)d‘^ ^^Whdx 


K&Th 


dx'^ 

K&Th'^ ^ ® l<j<daeMij 


X H4,K\d^-^^^h\o.K. 

KGTh 

where we have used the fact that maxmaxi3/^(j, a) < 

j x£K 

A further application of inverse estimate yields 

X < Ch"^ X \'^U,K\y^^hWK 

KdTh K( 3 Th 

< Ch?'\u\i^Q,\wh\2,h- 


(31) 


(32) 
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Then, we estimate the second term of 


L 


K 

i,2 








dxi 


^ ^ Bf{j,a){Ilo^K - Id)d°‘whdx 

j¥=i 


+ 



— Id)d°‘whdx. 


(33) 


According to dH]), since 



)^-/i|](no,if-/d)(ci+C2Xj) dx 


and 


0 , Cl, C 2 are constant coefficients, 



^j,c) 


3 


h'^ixj 


Xj,c)](no,ic - Id)c 3 dx 


0 , C 3 is a constant, 


thus, we can get that 



K — Id)d°‘whdx = 0. 


(34) 


The interpolation error estimate and inverse estimate yield 

f {Id-Uo,K)p- E Bf{j,a)iUo,K-Id)d^Whdx 

(35) 

< Ch^|u|4,_R-|V^W;i|o,if■ 


Then, we can get 

A combination of dsa and (EH) leads to 
E,Xu,Wh)= Y. 


Ken 


(36) 


(37) 


Similarly we obtain further, 

Eh{u,Wh) <Ch'^\u\4,Q\wh\2,h- (38) 

This, combined with the approximation error estimate, finishes the proof. □ 

Remark 3.2. Compared with and \1S^ . we prove the 0{h?) energy norm 
eonvergence rate without assuming the uniformity of the meshes. Besides, we 
only need the lowest regularity assumption u € i7^(f2). 
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3.3 Error analysis of the Adini element in norm 

In this section, we present the lower bound estimate of the error in norm. 
This is a generalisation of the result in [7] to arbitrary dimension. The main 
result of this section is the theorem below. 

Theorem 3.3. Letu and Uh he solutions of problem and i2S\) . respectively. 

Suppose that u € s > 4 and s > f + 1. Then, provided 

W/Wl^q) ^ 0, 

IIm - M;i||L2(n) >/?/^^ (39) 

where /3 = <5/||/||L2(n). 

Remark 3.4. By the embedding theorem of the Sobolev space, we need higher 
regularity of the solution in higher dimensions in order to guarantee H'^{K) C 
C^(K). Furthermore, it ensures the eontinuity of interpolation operators. 

Remark 3.5. For the rectangular domain it, the condition ||/||l 2 (q) ^ 0 implies 

I dfidxj IgMS'i) 7 ^ 0, 1 < i j < d. In fact, = 0, 1 < i j < 

d, then u is of the form 


d d 

M = ^ ^ CjjXjXj + y^^g{xi), 
i—1 i—1 

for some funetion g{xi) with respect to Xi. Then the boundary condition indi¬ 
cates u = 0, which contradicts u ^ 0. 

We postpone the proof of Theorem 13.31 after several technical lemmas. 
Define the global interpolation operator n?i to 14 by 

R/ilif = Uk for any K £%. (40) 

By means of Lemma 12.21 we can obtain the following crucial result. 

Lemma 3.6. For u G H§{n) P| H‘^{n), s > 4 and s > ^ + 1, it holds that, 

{Vl{u - Uhu), V^nhM)i2(n) > ah'^, (41) 

for some positive constant a, which is independent of the mesh size h provided 
that ||/||L 2 (f 2 ) 0 and that the mesh size is small enough. 

Proof. Given any element K, we follow [7] to define Pkv G Pi{K) by 


/ y^PKvdiX= [ V^vdx, l = 0,l,2,3,A, 
JK J K 


(42) 


for any v G i7®(D), (s > 4 and s > | + 1). Note that the operator Pk is 
well-defined. The interpolation operator Pk has the following error estimates: 


|r! - Pkv\j,k < Ch'^ ^\v\a,k, j = 0,1, 2, 3,4, 
\v - Pkv\j,k < Ch\v\j+i^K, j = 0,1, 2, 3, 


(43) 
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provided that v G H^{n), (s > 4 and s > | + 1). It follows from the definition 
of Pk in that 

V'^Pkv = Uo^K^^v. (44) 

By the aid of Pk, we have the following decomposition 

- n/iu), V^n?iU)i2(Q) = {VI{Pku -IlKPKu),VlllKu)i^2f^K) 

K^Th 

+ X! ~^K){Id - PK)u,V\liKU)L^(^K) 

K&Th 

= Ji + J2- (45) 

We first analyze the first term Ji on the right-hand side of (H5l) . By means 
of Lemma 12.21 the first term Ji on the right-hand side of (H5l) can be rewritten 
as 


j _ f O'^Pku O^Uku 


K^Th 

d 


Ik 


dxldx'i dx^ 


dx 


^ 3 


-- V V V ^ ^ ^ dx 

Q / a 2 2a 2'^-^ 
K&n i=l l<J<d dK OX^Xj OX, 


d 


EE E ? 


h] f d'^ild- PK)ud^IlKU 


K^Th 


Ik 


dx^dx^, dx^ 


dx 


EE E 

K&Th i=l l<i<d 


r d^u d^ild-UKu) 
y JkO^ 


dx. 


(46) 


Since |^L, = 0, 


du 


= 0. and Tpf^L, = 0, 


a-^u 


d'^xidxj . ’ d'^Xidxj If" 

integrating by parts yields 


= 0 , 


EE E 

K&Th j=l l<3<d 


hZ f d^u dh 


Ik dxfx] dxf 


= - E E E 


K^Th i=l l<i<d 

jZi 


d^u .2, 


IK 


= -EE E ? 


hZ.. d^u 1,2 


3 " dxfxj 

i—1 l<i<d * 

3¥=i 


By the commuting property of (14411 . 
d^{Id-PK)u 


dx^dx^ 


= {Id — no,x) 


d^u 

9x25x2 


, I <i j <d. 


Note that 


^i|92nKW|| I 

Z^ll Q^2 llL2(if) - 


i=l 
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This and the error estimate of (H51) yield 


K^Th j=l l<3<d 


(47) 

We turn to the second term J 2 on the right-hand side of (HHI) . By the Poincare 
inequality, and the commuting property of (l44ll . 


IJ2I < Ch^ Y. W^iild- PK)u\\L^^K)\uU,n 

K&Th 

< Ch^\\{Id-UQ)Vtu\\mn)\u\3,n. (48) 

Since the piecewise constant functions are dense in the space 

||(/d — no)V^M||i 2 (f 2 ) —>■ 0, when h —>• 0. (49) 

Summation of (ITTl) . (H51) and completes the proof. □ 

Again, the lemma below can be found in [7]. 

Lemma 3.7. Let u and Uh be solutions of problem h2U\) and \2S\) . respectively. 
Then, 


{-f,u- M/i)L=(n) = ah{u,nhu - Uh) - {f,UhU - Uh)L^{n) 

+ ah{u - UhU, u - Iffiu) + ah{u- lihU, Uh - Uhu) (50) 
-I- 2(/, UhU - m)l2(q) -I- 2ahiu - Utu, Utu). 


Proof of Theorem 13.31 It follows from (IHRl) that 

ah{u,IlhU-Uh)-if,TlhU-Uh)L^(n) < Ch'^\u\4^Q\IlhU-Uh\h < Ch^\u\lQ. (51) 
By the Cauchy-Schwarz inequality and the error estimate (I22L it yields 
ah{u-IlhU,u-Ilhu) + 2{f,IlhU-u)L2{a) < Ch‘^{\u\4^Q + \\f\\L2(Q))\u\4,n, (52) 

Ohiu - UhU, Uh - Uhu) < Ch^\u\l Q. (53) 

The error estimate of the last term of dsni) by Lemma 13.61 gives 

ah^ < ah{u — UhU, B/itt). (54) 


Hence, a combination of (ISmi - dSTl) leads to 

i-f,u-Uh)L^n)>Sh‘^. (55) 


for some positive constant 6, which is independent of the mesh size h and the 
mesh size is small enough. 

Therefore, 


\\u - M/i||L 2 (n) 


sup 

o/ii;ei2(n) 


{w,U- Uh)L^Q) 


> 


i-f,U- Uh)L^(Q.) 


- /IU=(n) 


>dl\\f\\L^(n)h^ 


This finishes the proof. 


□ 
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4 Numerical examples 

In this section, we present some numerical results of the three-dimensional 
Adini element by congruence partition of cubic meshes and non-congruence 
partition of domain to demonstrate our theoretical results. Herein, we give 
ui{x,y,z) = sin^{'!Tx)sin'^{Try)sin'^{Trz) and U 2 {x,y,z) = x^{l — x)^y^(l — 
y)^z^{l — as the exact solution of problem (l20l) . respectively. One can see 
the errors and the rate of convergence computed by uniform cubic meshes with 
the meshsize h = ^ tor some integer N in Figure 2. One can also see the errors 
and the rate of convergence computed by non-congruence meshes in Figure 3 
for logarithmic plot. 


The errors of u,(x,y,z) The errors of u^(x,y,z) 




Figure 2: The errors in and norms for ui{x,y,z) and U 2 {x,y,z) by 
uniform cubic meshes 


5 Concluding remarks 

In this paper, we studied the accuracy of the Adini element as a discretiza¬ 
tion scheme for fourth order problem in any dimension. We showed that the 
convergence rate is of 0{h‘^) order in energy norm in any dimension, and more¬ 
over, we show that the convergence rate can not be non-trivially higher than 
0{h?) order in integral norm in any dimension. By the Poincare inequality, 
we arrive at the conclusion that the convergence rate of Adini element for dis- 
cretising fourth order problem is 0{h^) in and the energy norm. This 

presents a complete exploration of the capacity of the scheme. 

The results provided in this paper are optimal in two-folded. On one hand, 
the full convergence rate of the energy norm is established under the assumption 
u € which is standard and of the lowest regularity. On the other hand, 

combining the two results together illustrates that neither of these two can be 
improved. The analysis of this paper has been sharp and economic enough. 
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The errors of u^(x,y,z) 



The errors of U2{x,y,z) 
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Figure 3: The errors in and norms for ui{x,y,z) and U 2 {x,y,z) by non¬ 
congruence meshes 


It is somehow surprising to observe the convergence rates in different norms 
are of the same order. This is because the Adini element function is not moment 
continuous across the edges, but it possesses internal symmetry on every cell 
when the grid is of tensor type. This point of view will hint us on other elements. 
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